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We develop a theory of magnetoresistance based on variable-range hopping. An exponentially 
large, low-field and necessarily positive magnetoresistance effect is predicted in the presence of 
Hubbard interaction and spin-dynamics under certain conditions. The theory was developed with the 
recently discovered organic magnetoresistance in mind. To account for the experimental observation 
that the organic magnetoresistance effect can also be negative, we tentatively amend the theory with 
a mechanism of bipolaron formation. 
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I. INTRODUCTION 

A. Organic magnetoresistance 

Organic magnetoresistance (OMAR) is a recently dis- 
covered large, low field, room-temperature magnetore- 
sistive effect {AR{B)/R < 10% a.t B = lOmT) in or- 
ganic light-emitting diodes (OLEDs) 0, |3|. To the 
best of our knowledge, the mechanism causing OMAR 
is currently not firmly established, although Prigodin et 
al. have put forward a model of OMAR based on 
spin-dependent carrier recombination and spin-dynamics 
caused by hyperfine interaction. For the benefit of the 
reader, we will briefly summarize some of the main ex- 
perimental results: 

1. It is an effect associated with the bulk resistance of 
the organic semiconductor layer 0,0 

2. It is only weakly dependent on the minority car- 
rier density , and occurs also in heavily p-doped 
devices 

3. It is independent of the magnetic field direction 

4. It obeys the empirical laws AR{B)/R oc B'^/{B'^ + 
Bl) or AR{B)/R (x B^ /{\B\ + Bq)^ dependent on 
material, where Bq ss 5mT in most materials 6] 

5. The Bq value is much larger in materials with 
strong spin-orbit coupling |3| 

6. It is only weakly temperature dependent 

7. It typically decreases with increasing voltage and 
carrier density 0,0 

8. Its sign can be positive or negative, dependent on 
material and/or operating conditions of the devices 
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Observation ^ establishes that OMAR is not related 
to carrier recombination in distinction to the model put 
forward by Prigodin. |SJ| shows that OMAR is caused by 
spin-dynamics. 



B. Correlation effects on hopping conductivity 

Although a solid understanding of carrier transport in 
disordered organic films has not yet been achieved, it is 
commonly believed that it occurs through variable-range 
hopping 00. Furthermore, it is known that Coulomb- 
correlation effects are much more important in organ- 
ics than in inorganics because of the smaller dielectric 
constant as well as confinement of the wavefunction to 
individual molecules. Electron-electron interaction ef- 
fects on hopping conduction have been studied for many 
years. For example, the intersite Coulomb interaction 
fiol IT]| results in the concept of " Coulomb gap" . Ya- 
maguchi, Aoki and Kamimura studied the effect of in- 
trastate Coulomb interaction (Hubbard-like interaction) 
in the strongly localized regime [l^ . Their model is based 
on the following Hamiltonian: 

II ^ ^ ^aa'^aa ^ ^ ^ U ct'^aa'^a — cr (1) 

Cqct, Ua and Uaa represent the one-electron energy, the 
Hubbard interaction energy and a number operator of 
state a with spin cr, respectively. Their model assumes 
a constant value for U and a uniform distribution for 
tan with a density of states v. The conduction mech- 
anism is assumed to be a one-electron hop assisted by 
one phonon. In the presence of the intrastate interac- 
tion there exist three electronic states for each Ander- 
son localized state, namely the unoccupied (UO) state, 
the singly occupied (SO) state and the doubly occupied 
(DO) state. SO states are found just below the Fermi 
level, Ep^ down to Ep — U and carry free spins, whereas 
DO states involved in hopping, which are spin-singlets, 
have an occupied deep level located near {Ep — U) and 
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accommodate an additional electron of opposite spin in 
a state located close to Ep (see Fig. QJ. Accordingly, 
there exist the following four different kinds of hopping 
processes: (1) hopping from SO to UO (2) hopping from 
SO to SO (3) hopping from DO to UO and (4) hopping 
from DO to SO, where the occupancy refers to that prior 
to the hop. In our treatment in section^we will neglect 
hops out of DO sites because of their smaller density of 
states (see Fig. QJ. YAK showed that Mott's law for the 
temperature dependent resistivity, p(T) still holds even 
in the presence of intrastate correlation : 



P{T) 
To 



Pooexp 
13.6 



(2) 
(3) 



where ^ is the localization length and k is Boltzmann's 
constant, although the value of ^ will differ whether U is 
considered or not. We have written Mott's law for the 
three-dimensional case, which may not apply to organic 
semiconductors. The replacement 1/4 ^ l/{d+l) in the 
exponent must be made in case of a different dimension- 
ality, d. Eq. 12 is valid only in the low field regime. In 
OLEDs, however, the apphed electric field is quite large, 
typically F w lO^V/ cm. For large fields eq. |21 remains 
valid, if the replacement T — ^eF/k is made, where e is 
the elementary charge [l3|. 
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FIG. 1: Schematic picture of the impurity states (left panel) 
and their density of states (right panel). Hops are allowed 
(solid arrows) from SO to all UO sites within the thermal 
energy kT from the Fermi energy, Ep = 0, hops between SO 
sites are allowed only for antiparallel spins. However, in the 
case of rapid spin-fiips hops between other SO sites become 
allowed (dashed arrows). Ec is an energy characterizing the 
density of states distribution below Ep- 



C. Magnetoresistance 

Although the intrastate interaction does not change 
the form of p(T), Kamimura and Kurobe showed that U 
plays an essential role in magnetoresistance When 
a magnetic field is applied, the spins of SO states be- 
come parallel to the field. Therefore, the probability of 
finding a pair of antiparallel spins becomes smaller with 
increasing magnetic field, and the hopping processes from 
SO to SO states become suppressed. Such a strong pos- 
itive magnetoresistance, which saturates at high fields, 
has indeed been observed experimentally [isl (l6l| . It 
has also been measured in strongly localized magnetic 
two-dimensional electron gases and interpreted using a 
similar model jl7|. However, the characteristic magnetic 
field scale above which the effect saturates is very large, 
orders of magnitude larger than the 5mT characteristic 
magnetic field scale of OMAR. Therefore, in this form 
the Kamimura-Kurobe model cannot account for OMAR. 
However, this model does not take into account spin-fiips. 
In organic semiconductors, it is well-known that there 
exist spin-flip mechanisms, the best-known of which is 
hyperfine interaction p^ . Hyperfine interaction results 
in the precession of the electron spin around the nuclear 
spins of the hydrogen atoms, and the precession period 
is typically 1-lOns. Furthermore, since the drift mobility, 
/i in disordered organic semiconductors is typically very 



low, about IQ^^cm^iys)^^ , the distance covered by the 
carrier within the precession period of the spin is only 
0.1 to Inm y (see note 0|), similar to the molecular di- 
mensions. Therefore, it should be necessary to take the 
spin-precession into account. In the limit of rapid spin- 
precession, hops become allowed to all SO sites for all 
carriers, rather than only for ones in a singlet state. Fur- 
thermore, since the hyperfine interaction is very weak, 
upon application of a magnetic field of typically several 
milliTesla the spin-dynamics is suppressed Q . 

The problem of variable-range hopping for large U, fur- 
ther taking account of the possibility of spin-flips, was 
considered by Meir and Altshuler 20] . The effect of spin- 
orbit coupling on hopping magnetoconductivity was also 
investigated by Shapir and Ovadyahu, where it was con- 
cluded that backscattering known from weak-localization 
theory is still important in the strongly localized regime 
[2ll |. Meir and Altshuler's calculation was only performed 
for a large applied magnetic field, which polarizes all SO 
states, thus blocking all SO to SO hopping processes. 
This leads to an effectively reduced density of states, and 
an exponentially enhanced magnetoresistance. There- 
fore the only possibility of such hops is through a spin- 
ffip with probability Psf- Their main result is that the 
variable-range hopping resistance is given by the parallel 
circuit of two resistors, pi and p2 given through: 



II. THEORY OF ORGANIC 
MAGNETORESISTANCE 

Sheng et al. 4] have shown that in organics 



PsFiB) = Psf(O) 



Br 



(6) 



where Bq is the hyperfine coupUng strength, or more gen- 
erally the strength of the spin-flip mechanism. Therefore 
application of B reduces Psf, because the associated Zee- 
man energy pins the spins, and it is straightforward to 
show that the theory by Meir and Altshuler can only re- 
sult in positive magnetoresistance. This is contrary to 
experimental observations for OMAR, which can be pos- 
itive or negative. However, we will pursue this theory 
further. At the end of our treatment, we will suggest 
a modification to the theory such that it can result in 
negative magnetoresistance also. 

The basic idea underlying the theory of variable-range 
hopping is that the average hopping distance is directly 
related to the density of states of potential target sites. 
Therefore, we need to determine the density of states 
available to a hop out of a particular site (site 2 in 
Fig. The density of UO target sites is determined 
by the density of states around the Fermi energy, vlEp), 
while that of SO target sites is determined by ^{Ep — U). 
Without spin-flips the overall density of target sites is 
V = v{Ef) + 1/2v{Ef - U) since half the SO states will 
have a spin opposite to that of the hopping carrier. With 
spin-flips however we have v — v{Ef)-\-v{Ef~U). Quite 
generally we can state that v — v{U,Psf)- Therefore 
we assume that Mott's law still holds, but the replace- 
ment V 1^(11, Psf) must be made in eq. 13 Finally, 
we note that a similar mechanism should also hold for 
nearest neighbor hopping, because of the possibility that 
the nearest neighbor site is an SO site with parallel spin. 
Then hopping into this site would be forbidden unless 
rapid spin- flips occur. 



III. DISCUSSION 

To simplify further discussion, we linearize the modi- 
fied Mott's law to obtain for the magnetoresistance: 



Ai?(g) 
R 

Furthermore, 



To\~ 1 Ai^(B) 
Y] d+l V '■ 



(7) 



where, in the second equation we have made the as- 
sumption, as it is commonly done in disordered organics, 
that v{E) = exp{{E-EF)/Ec) ^ (Fig.Cl right panel), 
although a Gaussian density of states is commonly be- 
lieved to be more accurate. If we moreover employ a 
crude estimate for Tq/T « (Ec/kT), then we obtain the 
following approximate result, which is useful for quanti- 
tative estimates and for examining qualitative trends: 



AR{B) 1 



R 



Q.d+1 



2{d+l) 



'Psf{0) 



B^ 



B^ 



Br 



(10) 



where a — Ec /kT . We would like to remind the reader 
that a is also the exponent of the current-voltage (I-V) 
characteristics predicted by the theory of space-charge 
limited current in the presence of traps, i.e. / cx 
which is routinely applied to OLEDs 22] . 

Unfortunately, we are not aware of any experimental 
measurement of U . It is important not to confuse U 
with Ui ~ beV used routinely in quantum chemical cal- 
culations for TT-conjugated molecules 0|, which refers to 
the interaction between electrons located on the same in- 
tramolecular site. However, we can use J7i to estimate 
U ~ ?7i(a/^)^, where awl — is the size of an in- 
tramolecular site. We therefore expect that U is about 
two orders of magnitude smaller than Ui. Finally, we 
estimate Psf(O), which will be determined by the com- 
petition between spin-flips and hops, therefore 



P5f(0) 



g^iBBo/h 



^iF/C + gfiBBa/h 



(11) 



A. Predictions of the model 

• Since the model is based on hopping conduction, 
OMAR should not exist in organic single-crystal 
devices with band transport. To the best of our 
knowledge, this has not yet been tested 

• The magnitude of OMAR depends critically on the 
ratio U/ Be- 
rn Our model is obviously consistent with experimen- 
tal observations (jTJ, ©, O, ^ 

• Our model is an example of a spin-dependent 
process which does not require thermal spin- 
polarization, which is in agreement with © 
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• In the high-field case the effect is predicted to de- 
crease with increasing F, since F replaces T in 
eq. Furthermore, F reduces the effect through 
eq. These statements are in agreement with O 

• The effect tends to be smaller in high-mobility ma- 
terials @. Theoretically this is a consequence of 
eq. ^2 and because there should also exist an in- 
verse relationship between fi and Ec- 

• Making a guess for U O.leV, taking typical val- 
ues for Ec « O.leV Q and a « 10 JJ, assuming 
rapid spin-flips Psf = 1 and d = 3, we obtain 
as an estimate AR/R « 10%, certainly consistent 
with OMAR 



B. Speculative modification to achieve negative 
magnetoresistance 

The main shortcoming of the model is that it predicts 
only positive magnetoresistance, in contradiction with 
experimental observation |(SJ). Therefore, a modification 
of the theory is required. Although we have been unable 
to find a mechanism that is firmly established, we are able 
to suggest a candidate mechanism for which some experi- 
mental evidence exists. As we described above, spin-flips 
will lead to an increase in the population of DO (spin 0) 
sites at the expense of SO (spin 1/2) sites. We propose 
that bipolaron formation [23 can take place at DO sites. 
The formation of bipolarons is well-established for heav- 
ily doped organic semiconductors H^HJ^^, evidence in 
OLEDs is much sketchier ■ Bipolarons are com- 

monly expected to have lower mobility than polarons, 
which are well-known to be the spin-1/2 charge carriers 
in organic semiconductors. This expectation results from 



the observation that, to account for the binding energy 
of the bipolarons, the lattice distortion around a bipo- 
laron is greater than that around a polaron, resulting in 
additional mass and reduced hopping probability. Since 
the applied magnetic field reduces spin-flips and therefore 
occupation of DO sites, bipolaron formation will be sup- 
pressed, resulting in negative magnetoresistance. More 
generally, any mechanism through which the correlation 
length of DO sites becomes less than that of SO sites can 
result in negative magnetoresistance. 



IV. CONCLUSION 

We have proposed a theory for the recently discovered 
organic magnetoresistance effect. Our theory is based on 
the variable range hopping model which states that the 
average hopping distance is determined by the density 
of states of potential target sites. We show that, in the 
presence of Hubbard interaction, this density is restricted 
by Pauli's principle unless rapid spin-flips occur. In or- 
ganics, spin-flips are known to result from the hyperfine 
interaction, which can be eliminated by the application of 
a weak external field. Therefore the application of an ex- 
ternal field leads to an exponential increase in resistance. 
It has been brought to our attention that Bobbert and 
coworkers are currently developing a similar model. 
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